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ABSTRACT 

In the last few years before merger, supermassive black hole binaries will rapidly inspiral and 
precess in a magnetic field imposed by a surrounding circumbinary disk. Multiple simulations suggest 
this relative motion will convert some of the local energy to a Foynting-dominated outflow, with a 
luminosity ~ 10^^ ergs^^ (i3/lO*G)^(M/lO^M0)^(?;/O.4c)^, some of which may emerge as synchrotron 
emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. 
On top of a secular increase in power (and v) on the gravitational wave inspiral timescale, orbital 
motion will produce significant, detectable modulations, both on orbital periods and (if black hole spins 
are not aligned with the binary's total angular momenta) spin-orbit precession timescalcs. Because 
the gravitational wave merger time increases rapidly with separation, we find vast numbers of these 
transients are ubiquitously predicted, unless explicitly ruled out (by low efficiency e) or obscured (by 
accretion geometry fgeo)- If the fraction of Poynting flux converted to radio emission times the fraction 
of lines of sight accessible fgeo is sufficiently large {fgeo^ > 2 x 10"'' for a 1 year orbital period), at 
least one event is accessible to future blind surveys at a nominal 10^ deg^ with 0.5 mJy sensitivity. Our 
procedure generalizes to other fiux-limited surveys designed to investigate EM signatures associated 
with many modulations produced by merging SMBH binaries. 

Subject headings: black hole physics — cosmology: observations — radio continuum: general — surveys 



1. INTRODUCTION 

Merging supermassive black hole binaries should nat- 
urally possess a circumbinary accretion disk whose MRI- 
driven turbulence generates and imposes a substantial 
external magnetic field. Recent simulations suggest 
that supermassive black holes (SMBHs) moving through 
an imposed magnetic field should pro duce a Foynting- 
domina ted outflow, consisti ng of a jet iFalenzuela et al.l 
(|2010a[) and diffuse emission iMoesta et all (|2011[) . Thu"s? 
even in the absence of accretion, during their last 
years before merger, SMBH binaries should generally 
have faint emission, modulated by their orbital motion. 
Though the outflow (henceforth "jet") power increases 
d uring the inspiral, end ing in a bright flare, as discussed 
in I Kaplan et al.l (I2011D , such flares are too faint or too 
rare to easily detect, unless almost all the Foynting flux 
is efficiently converted to low-frequency radiation. On 
the other hand, because SMBH binaries spiral in very 
slowly through gravitational wave emission, each merger 
fiare is preceded by a long phase during which the jet 
is modulated and growing, at only slightly reduced effi- 
ciency. These modulations should be easily accessible to 
future radio surveys. 
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In this paper we calculate how frequently modula- 
tions will occur and how often they can be detected. 
O ur paper adopts and extends the assumptions used 
in iKaplan et al.l ()2011[ ) for the jet power versus binary 
masses and for fiducial SMBH binary merger rates. In [12] 
we sum over all SMBH binaries, to determine the num- 
ber of binaries on our past light cone with detectable 
modulation. In ^we describe how long, how frequent, 
and how bright modulations from each SMBH binary 
should be. Finally, 311 we discuss, in order to detect the 
signature presented here, targeted EM and GW surveys 
must overcome limitations of intrinsic AGN variability 
or small limiting distances, respectively. We also briefly 
explain how the many distinctive variations in the pre- 
dicted light curve will distinguish this source from other 
candidate modulation. 

To date, most 1 GHz radio transient surveys have sur- 
veyed at most a few tho usand square deg rees with rms 
sensitivities of ~ rnJy (jOfek et al.l l201ll l. with many 
having smaller area, less sensitivity, and sporadic sam- 
pling. Upcoming surveys will be larger, more regular, 
and more sensitive. I n particular, the VAS T (Variables 
and Slow Transients; iMurphv et al.l ()201lD ) project us- 
ing the A ustralian Square Kilom eter Array Fathflnder 
(ASKAF; iJohnston et all (l2007h 'l wiU survey roughly 
10^ deg^ daily down to a nominal sensitivity of 0.5 mJy. 
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We adopt these parameters to motivate our discussion, 
keeping in mind that survey plans often change and 
that unknown factors such as efficiencies or obscura- 
tions strongly impact our results. Rather than fix these 
choices, we leave in scalings, so interested parties can 
make their own predictions. 

1.1. Context 

Several mechanisms for precursor, prompt, and 
delayed electromagnetic signatures of SMBH merger 
have been propose d , inc lu ding Poynting-do minated 
iets (IKaplan et al.l 120111: IPalenzuela et all I2010bl: 
iMoesta et al.l 1201 1[ ): disk emissi on powered by post- 
merger 



200,^ _ 
Rossi et alT 



perturbations a nd sh oc ks (ISchnittman fc KrolikI 
lO'Neill et al.l l200l IMeeevand et al.l 120091: 



120101) 



^ , t he viscous-driven refilling of 

the post-merger disk (jMilosavlievic fc PhinnevI l200l 
iTanaka fc Menoul 120101 ): accretion of a fos sil inner disk 
drive n by the smaller body's inspiral (iChang et al 



201C ) ; increased tidal di sr uption 



200a IWegg fc Bodl I20I0I: IStOT^ 



rates 
: Loebl 



( Chen et al 
2011D and 



perturbations of galactic-center stellar cores: and direct 
modul at ion of a clo s e circumbinary disk (Haiman et al.l 
2009bt lBode etal.1 [20TTI : iRoedig et all 1201 ll ): see 
SchnittmanI (I2010D and references therein. These 



merging binaries will also produce gravitational waves 
potentially a ccessib le to pulsar timing arrays and LISA 
ilJenet et al.l I2009t Ipemorest et all 120091: ISesana et all 
120091: IHugh cs"2002l: ISesana et al.ll2007[) . However mea- 
sured, the SMBH merger rate would provide invaluable 
direct constraints on the merger ra tes of galaxies (versus 
direc t methods, as discussed in iBertone fc Conselicd 
120091 and references therein) and the form ation and 
evolu tion of supermassive black holes (see, e.g. iVolonteril 
I2OIOI and references therein). Depending on the identi- 
fication epoch, mechanism, and followup observations, 
direct detection of merging SMBH binaries could pro- 
vide infomation about accretion, stella r dynamics, the 
evo lution of galaxies (jVolonteril 120101 : iBloom fc et all 
l2009i) . 

A number of binary SM B H are known or sus- 
pected (e.g . . iKorrioss al l2006t iRqdriguez et all l2006al : 
ISmith et all 120101: IBurke-Spolaod I2010D . While many 
dual AGN candidates (e.g., 3C75 and NGC 6240) 
have been discovered with kpc-scale-separation, very 
few candidate SMBH binaries are known with < 
10 pc separation. Objects with shifted broad line re- 
gions like SDSS .1092712.65-^2 94344.0 (IKomossa et al l 
I2OOI or J153636. 22-^044127.0 (iBoroson fc Laueiil2009D 
may be SMBH bin ar ies: s ee, e.g.. iBogdanovic et al.l 
(I2009D: iDotti et al.l (I2009D and references therein. 
iRodriguez et al.l (|2006b( ) demonstrate 0402-^379 has 
two resolved cores with a projected 7pc separa- 
tion. More speculatively , OJ287 has precessing jets 
(jMarscher fc JorstadI 120111 ) and outbursts with an ap- 
proximately 12 year period ([Sillanpaa et al.ll988() . These 
modulations have b een fit both by binary motion 
(jValtonen et al.' '2008|) and by a warped accretion disk 
(see, e.g. iKatz, .1997i r Since most SMBH binaries are 
too far away to be resolved by VLBI, future candi- 
dates will likely be found by o ffset emission lines or 
modulated spectral features (e.g..lCoinerford et al.ll2009l : 
iShen fc Loebl 120 10'-: 'Tsalmantza et alj [20111). For exam- 
ple, recentlv. lEracleous ct al.. (,20 IL) found 14 quasars 



with offset H/3 emission lines, with offsets that shifted sig- 
nificantly between widely separated observations. These 
deviations could arise from line-of-sight orbital accelera- 
tion in an SMBH binary. 

The search for binary SMBH is complicated by AGN 
variability. This variability is largely not periodic; de- 
tecting such periodic behavior in a radio light curve 
would be a strong indica tion of an inspiralling SMBH 
pair (e.g.. lKomossall2006f ) . 

Our study complements the previous model-neutral in- 
vestigation by [Haiman ct al. (2009a), who also both cal- 
culate the number of SMBH binaries per unit orbital 
period and the conditions under which (optical) surveys 
could detect their variability; see, e.g., their Figure 9. 
They examined both gravitational-wave-dominated and 
disk-driven inspiral. By contrast, in this paper we study 
the gravitational-wave-dominated phase only; adopt a 
concrete emission model; employ a state-of-the art pop- 
ulation of merging SMBH binaries distributed over cos- 
mic time; and address how AGN background variability 
and jet beaming would impact survey performance. We 
also discuss how gravitational wave searches with pulsar 
timing naturally complement EM surveys for massive, 
comparable-mass SMBH binaries. 

2. POPULATION OF MODULATED OUTFLOWS FROM 
MERGING BINARIES 



Foll owing IPalenzuela et al.l (|2010b( ) and IKaplan et al.l 
(I20T1I) . we assume each merging binary black hole has a 
Poynting-dominated outfiow, even in the absence of ac- 
cretionQ For convenience, we scale the outfiow's instan- 
taneous luminosity to a fraction eEdd of the Eddington 
luminosity at the merger event: 



Lflare — ^Edtiq'^ {v / Vmax)'^ L, 



dd 



(1) 



where q = raijrax < 1 is the binary's mass ratio and 
w = /3c is the binar y's coordinate re l ative c ircular ve- 
locity {\dt(ri - rx)\. IPalenzuela et al.l (|2010bD find that 
most of the emission is trapped into collimated jets, with 
^Edd ~ 0.0 02. Using a r efined calculation of the out- 
going flux, IMoesta et al.l (|2011| ) calculate that the jets, 
while present, are associated with a much brighter dif - 
fuse quadrupolar emission; in total, IMoesta et all (|201lD 
estimate the same fields correspond to eEdd ~ 0.02. In 
either scenario, the Poynting-dominated outflow will nat- 
urally interact with the strong ambient magnetic field 

^ On physical grounds we expect Ljet oc M^q^v^ , based on 
the black hole area, the energy density of the magnetic field, and 
the relative orbital velocity of the hole and ambient field. Self- 
consistent MRI amplific ation already prod uces a large enough B 
field in the disk (see, e.g. IPessah et al] l2006') to p ower the instanta - 
neous outflow assume d abov e. Following Bcck with et al.l II2009I ). 
IRothstein &: Lovelace! I|2008l ). and references therein, we antici- 
pate small loops will escape into the circumbinary region, expand 
rapidly in the absence of confining material, and form a fluctuating 
but partially coherent large-scale flcld in the interior that threads 
the circumbinary region and the hole. Ordered ambient flux can 
be advectod inward, increasing the ordered flcld near the horizon. 
Neither process requires matter accretion. Still higher magnetic 
fluxes can be achieved in the presence of accreting matter, a small 
amount of which will flow inward from the nonrelaxed inner edge of 
the circumbinary disk. A detailed discussion of B field evolution 
and flux advection onto the horizon is beyond the scope of this 
paper. Rat her than adop t a mo del for B holds near the hole, as 
described in lKaplan et al.l I I2011I) we conservatively adopt a smaller 
magnetic field, such that the maximum jet luminosity just prior to 
merger is the Eddington luminosity if eEdd = 1- 
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and coronal plasma to produce radiation. Lacking a de- 
tailed emission model, we assume a fraction eiadio of the 
beam power is converted to radio as F,, cx emdioL/iy 
and absorb all details of spectral model, K-corrcction, 
and most issues pertaining to anisotropic cmissior0 into 
Cradio- A survcy with a flux threshold Fn^rnin will there- 
fore be sensitive to all binaries inside a luminosity dis- 
tance dEdd('y, -M7 l)'- 



: 14.2 Gpc(/3//3„ 



/92(M/106Afo)(e/0.002) 



(^^.,mi„/niJy)(;^/GHz) 



: ^radio^Edd 



,(2) 
(3) 

where M is the total mass of the binary. All key re- 
sults and figures are presented as functions of e, to al- 
low efhcient scaling to any physically motivated values 

of eEdd, Eiadio- 

The jet power gradually increases over time as the bi- 
nary spirals inward towards merger due to the influence 
of gravitational radiation, with circular velocity increas- 
ing as 



(4) 



dt Gihivr 

where the symmetric mass ratio 77 = mim2/M'^ = q/(l + 
q)^. As each outflow is tied to a single black hole, we 
anticipate modulation at the orbital period 



C 

P=—2tiMP 



-3 



(5) 



Finally, for black holes with spin S = ^to^ misaligned 
with the orbital angular momentum, the orbital plane 
should precess, leading to mod ulations on a timescale of 
order (jApostolatos et aLllTOOl ) 



27r 



G 



|^(2 + 3g)c.2M^6 
^(2 + 3^)^2^/35 



(late) 
(early) 



(6) 



In these expressions wc differentiate between early in the 
inspiral, when the orbital angular momentum is much 
larger than the black hole spin, and late in the inspiral, 
when L < Si. Though the latter case occurs only in the 
last few orbits before merger with significant spin-orbit 
coupling, it ensures a significant change in the direction 
of L during each precession period and provides the best 
opportunity for jet modulation. On the contrary, early 
in the evolution, the orbital plane precesses through a 
small angle ~ 0{S/L); while many cycles occur, only 
for exceptionally bright jets with well-understood vari- 
ability could we ascribe small modulations to precession. 
Spin-orbit misalignment may not occur in gas-rich merg- 
ers, as accretion may a lign spins to the orbital plane 
(|Bogdanovic et al.ll2007f ). 

Given these timescales and a birthrate dtn per unit co- 
moving volume and time, we can determine the number 

^ Later we will explicitly account for how orientation-dependent 
obscuration (fgeo) or beaming {/beam) impacts the number versus 
flux distribution. 



of SMBH binaries on our past light cone whose modula- 
tion timescales lie in a desired range and whose jets are 
bright enough to detect: 



Ni = I dV, 

dVcdtn [T(/3+ 



dp 

f}^(z) df3/dt 



TiP-)] 



(7) 



where dVc is the comoving volume element, where T{v) 
is the time until merger for a binary of velocity c/3 



cy 2,0b 1] 



(8) 



and where the redshift-dependent expressions v± are 
defined using the smallest interval satisfying (a) 
the desired time-velocity relationship; (b) the con- 
straint that d < d^dd', (c) that T{f3±) are less 
than the time since decoupling from the circumbi- 
nary disk - 0(1) Myrq7A3(jVf/io8jy^^^i7/i3 (gg^^ ^ 

ISchnittman fc fcolikl 120081 : lO'NeiU et al][2?M and ref- 
erences therein) For example, the total number of bina- 
ries whose observed periods lie in [Pobs,~i Pobs.+] follows 
from 



I3c = I3,nax \l 47rd|F™i„ Aiz/Ledde?^ (9a) 

/3H. = r^i(0.86Myr x q'''^^{M /IQ^M^f'^^) (9b) 

/3± = min[/?™,„ max[p-i(Pob,.^/(l + z)), /3„ /3,]](9c) 

where T^^,P^^ are the inverse functions of Eqs. I5I8I a 
similar expression with P — > Tproc holds for precession. 
This expression distinguishes between the source-frame 
period P and the observed period Pots = P(l + z)- Fig- 
ure [T] illustrates how these four critical surfaces define 
a range of orbital velocity /3 consistent with maximum 
source distance, maximum lifetime, and allowed period 
range. For example, circumbinary disks limit the orbital 
period of decoupled, GW-inspiral- dominated binaries to 
P < 3yrqi5/26(^,j/;^Q8^^g)29/26^ Similarly, circumbinary 
disks a priori limit the number of decoupled BH binaries 
to less than ~ 1 Myr times the all-sky merger rate: less 
than 10^, or < 240deg~^. This limit is reached only for 
very sensitive surveys that probe all black holes through- 
out the universe. As implied by Figure [TJ for typical sur- 
veys and for the vast majority of decoupled binaries at 
moderate redshift, the flux limit or period limit bounds 
the space of detectable binaries; the circumbinary disk 

plays no role; 

Following iKaplan et al.l (j2011| ). we calculate the distri- 
bution of modulated flux at the earth using a distribu- 
tion of merger rates as a function of black hole mass 
and cosmic time. The assembly of SMBHs is recon- 
structed through dedicated Monte Carlo merger simu- 
lations which are framed in the hierarchical structure 
formation paradigm. Briefly, these models evolve the 
BH population starting from BH "seeds," through ac- 
cretion episodes triggered by galaxy mergers, and in- 
clude the dynamical evolution of SMBH-SMBH binaries. 

^ For brevity, in our expressions we assume an accretion rate 
rh = O.SMgjj disk viscosity parameter a = 0.3, following 

IKing et al.l mOT\ ). 
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The SMBH population is consistent with observational 
constraints, e.g., the luminosity function of quasars at 
1 < z < 6. the M — a relation and the BH mass density at 
z = (jVolonteri et a l. 2003, 20081 iVolonteri fc Begelmanl 
l2010f l. To illustrate the situation, we ad opt two of 
the fid ucial merger distributions as used in lArun et al.l 
(|2009| ). The two models used here are representative of 
two proposed seeding scenarios (see iSesana et al.l 12011a 
for fu rther discussion). In the notation of lArun et al. 
(|2009h . we compare models LE and SE, where S versus 
L refers to the seed size - large or small - and E refers 
to "efficient" accretion. The models described in that 
paper are representative of a range of plausible SMBH 
growth scenarios. As with uncertainties in Eradio, we at- 
tempt to make our conclusions robust to specific merger 
assumptions. 

Any reasonably-large sky patch should contain many 
bright SMBH binary jets whose plausible modulation 
timescales are within reach: minutes to months. Pre- 
cessing sources, though the most unambiguous sources of 
modulation, are rare: only for the most optimistic effi- 
ciencies (e ~ 1) will a typical-scale survey (10** deg^) have 
sufficiently bright precessing jet in its sky patch [Figure 
[2]. On the contrary, many binaries have orbital periods 
between minutes and months and are close enough to 
produce bright jets [Figure [3] . The slow nature of gravi- 
tational wave losses insures most binaries are discovered 
near the widest orbits. By contrast, though the largest 
and closest binaries are most likely to be detected, any 
binary not inconsistent with the flux limit of the survey 
and that can plausibly be produced at that mass and 
redshift has a reasonable chance to be recovered [Figure 
|4]. On the one hand, each value of F/e defines a mini- 
mum mass versus redshift Mc{z) oc F/e below which no 
emission is visible even at merger|3 these limiting curves 
are shown in Figure |4] for several choices of F/e. On 
the other hand, at each redshift the merger process does 
not efficiently produce SMBHs above a critical mass. In 
our two merger trees, this maximum mass Mmgr{z) is 
approximately 

Mmar{z) ~ 1O1°-^/2-5M0 (10) 

As a concrete example, in Figure |4] we show the mass- 
redshift distribution expect ed for a fiducial surv ey with 
F = 0.5 mJy sensitivity (e.g. lMurphv et ani2011| ). adopt- 
ing the most optimistic efficiency e = 1. Though concen- 
trated at moderate mass and redshift, the distribution 
extends throughout the region bounded by the smallest 
fiux (black) and largest mass (black dotted) Results at 
other fiux sensitivities are qualitatively similar, derived 
from this figure by suppressing binaries near and beyond 
the desired fiux limit. Finally, the mass-redshift distri- 
bution implies a cumulative flux distribution A''(> F/e), 
with each value of A^(> F/e) corresponding to the inte- 
gral in mass and redshift above a specific threshold curve 
for F/e (blue). In particular, this relationship and these 

This mass-redshift boundary is defined by aLi^fifil A.'wd^ — ^mzn 
and would be attained only at merger of a comparable-mass binary. 

^ For numerical reasons, these boundaries also determine the 
range of fluxes we can reliably model. Our catalog consists of 
redshift bins Az ~ 0.1 out to z ~ 10. For fluxes F/e < 0.1 mJy, 
binaries from z > 10 may contribute significantly. By contrast, 
for F/e > 10"* Jy, the limiting sensitivity passes inside the smallest 
redshift bins for a potentially significant proportion of mergers. 



contours suggest that for every flux limit F/e < 10^ Jy, a 
significant fraction of all deteted binaries are always high 
mass and moderate- to high-redshift sources {z > 1). 

Relatively few AGN have unobscured lines of sight to 
their central engines, particularly at radio frequencies. 
Adopting a parameter /goo to characterize the fraction 
of unobscured lines of sight Q the cumulative distribu- 
tion of sources versus flux is fscoN{> F/e). Fitting to 
the brightest sources in Figure [3l we estimate the num- 
ber of detectable binaries with periods less than 1 yr is 
betweerQ 

A^obs^lO-30/gco (^^) (11) 

for a moderate range of flux {F/e G [10~^, 10^] Jy). Tak- 
ing the lower limit, a fiducial survey with 10^ deg^ cov- 
erage and 0.5 mJy sensitivity should have at least one 
binary with bright, detectable jets and observed period 
less than 1 yr if the source obscuration and composite 
radiative efficiency satisfy 

/geoc > 2 X 10-4(^^™„/0.5 mJy) (12) 

In terms of Figure [3l this constraint corresponds to the 
number of orders of magnitude between (a) the flux at 
which the sky coverage limit (dotted) and the blue cu- 
mulative intersect and and (b) the survey flux limit (just 
off scale to left). The number of sources in wide orbits is 
large; only extremely inefficient conversion of energy to 
radio can prevent them from being seen. 

Conversely, Figure [2] shows only a few binaries are 
both massive enough to produce a detectable radio jet 
by this mechanism and tight enough to undergo pre- 
cession on an accessible timescale. Empirically Figure 
[2] is approximated by N{> F) ~ 10"i-^(F/e Jy)-2/3. 
Therefore, even an equally wide (lO^deg^) but much 
deeper (0.05 mJy) survey will have at least one binary 
with bright, detectable precessing jets only if jet power 
is efficiently converted to radio energy: 

fH^e > 0.2(^:™„/0.05 mJy) (precession) 

For brevity and clarity in this discussion we assume 
the outfiow emits roughly isotropically in all directions. 
If instead the same amount of power is beamed exlusively 
along a fraction fbeam < 1 of all lines of sight, the number 
of detectable sources versus fiux Nteam can be calculated 
from the previous expression via correcting for the frac- 
tion of the sources pointing towards us and the increased 
fiux emitted along the remaining lines of sight: 

Nbeami> F/e) = fteamN llt^^ (13) 

Ignoring detectors' sensitivity, the fraction of jets point- 
ing towards us and in the circumbinary phase is less than 

® This parameter absorbs all larger-scale obscuration effects, 
such as any obscurin g torus or c louds that appear in AGN uni- 
ficatio n schemes; see I U rrv fc Pad ovani (1995), Lawr ence fc Elv3 
II2OIOI ) and references therein. Standard AGN unification schemes 
require the torus subtend 65°, corresponding to a factor fgco — 2—4 
in the radio l|Lawrence fc Elvi3l2010l 'l. 

At the highest flux, the local universe dominates. Based on 
Eq. [7] and a uniform local merger rate, we expect that N{> F/e) 
at large flux must scale as Af oc M^I^P^I^ /(F /efl"^ . 



the total number of binaries that have decoupled from 
their circumbinary disk on our past light cone (< 10^) 
times fbeam = Ar2/47r, the relative solid angle covered by 
a jet. By contrast, for experimentally accessible and ob- 
scrvationally interesting flux ranges, the number versus 
flux trend iV(> i^) oc 1/F [Eq. [II] and Figure H lead- 
ing to Nbeam independent of beaming^ While strictly 
true, beam modulation allows much longer orbits to be 
accessible and thereby can still significantly increase the 
number of accessible sources, as described below. 

3. JET MODULATION 

Throughout the inspiral, each black hole drives a 
Poynting-dominated outfiow, powered by the strong am- 
bient fiel d provided by the c ircumbinary disk. As de- 
scribed in lKaplan et al.l ()2011[ ) (cf. their Footnote 1), our 
assumptions correspond to assuming the disk provides a 
magnetic field such that the jet luminosity near merger 
Lj(/3 = Pmax) oc B'^M'^ is the efficiency times the edding- 
ton luminosity (cEddiedd); the corresponding constant B 

field is B ~ B^dd =i 6 X lQ'^e]l^^{M /KfMQy^'^G. A de- 
tailed treatment of the outflow kinematics and emission 
driven by this ideal-MHD j et is beyond the scope of this 
paper; see, e.g.. Section 5 o flLvutikovl (120111) for plaus ible 
emission microphysics and iPalenzuela et al.l (|2010b( ) for 
MHD. However, the large circumbinary field and coro- 
nal plasma provide a natural mechanism for dissipating 
the outflow, through strong synchroton losses, inverse 
compton scattering off of accelerated particles, an am- 
bient medium to shock, and naturally relat i vistic par- 
ticle velocity scalcsPlFollowing iKaplan et al.l (j2011| ). we 
assume that a fraction eiadio of the jet power is eventu- 
ally emitted at radio frequencies, associated with a jet 
tied to the spin axes of each black hole. 

As the black holes orbit and precess, the flux along 
our line of sight will be modulated by orbital effects 
e.g., bolometric Doppler boosting, as in lLoeb fc Gaudil 
2003| )): by precession of the spin axes changing the ori- 
entation of the jet; and other relativistic effects associ- 
ated with an ultracompact orbit. The observable mod- 
ulation depends sensitively on jet dynamics and emis- 
sion mechanisms, all beyond the scope of this paper. As 
an example with efficient conversion (e < 1) and sim- 
ple modulation, in the rest of this paper we optimisti- 
cally assume Foynting flux is efficiently re-radiated by a 
sufficiently dense corona of weakly relativistic electrons 
as synchrotron radiation, in a short optically thin low 
Lorentz factor "jet" . In this scenario, modulations at the 
orbit period arise from Doppler boosting during the orbit 
of each unresolved jet, with relative amplitude ~ P/ Pmax 
times factors of order unity depending on the line of sight, 
spectrum, and mass ratio. Modulations at the preces- 
sion period arises as the jet flow direction relative to the 
line of sight rotates, with line-of-sight flux fluctuations 
of order unity for a low-Lorentz factor jet. On the con- 

* For surveys limited to very low-redshift sources, N{> F) oc 
p— 3/2 rp^ ^Yie extent that all beamed sources are also nearby, 
beaming would indeed increase the detectable number: Nij^am c< 
,-1/2 

beam ' 

^ For example, the Poynting flux in this idcal-MHD jet {E.B = 
0) has a strong electric field \E\ ~ y^Lj/r^ ~ B^^^{M/r). Even 
in the absence of direct acceleration, this field produces strong 
transient drift currents v ~ c. 
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Fig. 1. — Critical orbital velocities versus mass: For an 

equal-mass binary of mass M at z = 0.6, a plot of orbital velocity 
versus mass. The two dotted upward-sloping lines are the velocities 
associated with 1 year (bottom) and 1 day (top) orbital periods; a 
one year survey with one day cadence will be roughly sensitive to 
this range of time periods. The downward-sloping blue line shows 
the minimum velocity /3c, above which the jet is bright enough to 
be seen at z = 0.6, assuming e = 0.002 and F^in = 1 mjy; this 
curve depends strongly on source distance, efficiency, and survey 
flux limit. Finally, the thick green line is /3», the velocity at which 
the binary separates from the circumbinary disk [Eq. [9]. The 
shaded region is bounded by /3± and corresponds to the range 
of velocities where a binary is both isolated and detectable as a 
periodic source. The time AT a binary spends evolving through 
this region is AT = T(/3+) - T(^_) for T given by Eq. |S1 
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Fig. 2. — Sources with precession modulation: Total num- 
ber: Assuming all binaries have significant spin-orbit misalign- 
ment, the number of sources with precession periods Pmin = 1 vciui 
and Pmax = 1 yr (blue) or 1 month (green) over the entire sky, with 
flux above F/t. To guid e the eye to scales comparable to future 
radio surveys like VAST IIMurphv et al.ir2011l ) , a horizontal dashed 
line corresponds to 1/10* deg^; a vertical dashed line corresponds 
to 1 mJy. 



trary, models with highly relativistic outflows have much 
stronger beaming, much less efficient conversion of jet 
to radio power, and a weak if any tie between the radio 
modulation timcscalc and any binary period: few sources 
will point towards us, while each will be radio-fainter and 
more randomly variable. 

To identify modulation, the absolute change in flux 
from the modulation [F y. 5 for 5 the relative intensity 
change) should be larger than the detector noise. Equiv- 
alently, the per- measurement flgure of merit derived from 
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Fig. 3. — Sources with orbit period modulation I: Total 
number: At any instant of time, the number of decoupled SMBH 
binaries versus average jet flux, limiting to orbit periods within 
Pmin = 1 inin and Pmax = 1 yr (blue and dotted blue) or 1 month 
(green and dotted green) over the entire sky, with flux above F/t 
for e = Eradio^Edd the Composite efficiency of the jet, using the 
SE (small seed, efficient accretion; dashed lines) and LE (larg e 
seed, efficient accretion; solid lines) model of lSesana et al.l l l2011al ). 
For example, for an efficiency t = fO""^ and target flux sensitivity 
f mjy, one to several hundred binaries are bright enough to see 
and have periods less than one year; of order one has a period 
less than 1 month. In the limit of infinite sensitivity, all rates 
must be less than the all-sky rate (~ 30/ yr) times the longest time 
allowed (1 Myr, corresponding to the time since separation from a 
circumbinary disk). In the limit of poor sensitivity, the number 
versus flux scales as as usual for the nearby universe. For 

comparison, the solid black line is 10'^{F / ejy)~^f'^ . The top axis 
shows ticks at Zmax, the redshift at which L^^^/imij^ = F/t for 
M = W"^ Mq; as we assume e < 1, no source can be detected at 
higher redshift. 



Fourier transforming the recieved flux 
SNR = d-^ 



F„ 



(14) 



should be greater than 1. Though realistic surveys wiU 
have a higher detection threshold, this factor can also 
be absorbed into e. For Doppler boosting at the orbit 
period, the relative change in power is smah for bright 
and long-hved jets {q ~ 1), owing to minimal contrast 
between the two holes' jets: to order of magnitude the 
relative power increase S is [cf. Eq. [1] 



P 1 



(15) 



In this case, the analysis of 
with 



Pmax 1 ^" 

foUows, replacing Eq. ^ 



Pc,rel — Pn 



Ledd^q 



1 



1/3 



(16) 



Out of the population of binaries whose jet power is 
bright enough to be seen and an orbit period in a testable 
range [Figure [3], Figure [5] shows a not insignificant frac- 
tion have jet power bright enough that Doppler boosting 
could be accessible. By contrast, for precession-induced 
changes in the jet direction, we expect a sinusoidal rela- 
tive flux change of order unity: any precessing jet that's 
bright enough to be seen will have detectable modula- 
tion. 




Fig. 4. — Sources with orbit period modulation lb: Mass- 
redshift distribution: Contours of the mass-redshift distribution 
(solid black) of all binary black hole outflows with F/e >= 0.5 mJy 
(thick black curve) and observed periods Pobs between 1 yr and 
1 min. This scaled flux limit corresponds to the flducial VAST sur- 
vey sensitivity {F = 0.5 mJy) for the most optimistic conversion of 
outflow to radio power (e = 1). Contours are shown at 1/4, 1/2, and 
3/4 of the maximum value {dN/dlog A'Idz ~ 0.06). The distribu- 
tion shown is smoothed, built by convolving a gaussian kernel with 
the underlying merger tree, with smoothing lengths A log M ~ 0.3 
in mass and (A^: = 0.3) in redshift. Also shown are contours of 
the largest possible (Eddington-limited) flux at a given redshift: 
F/e = Ledd/47rd2 1/ for 10"^, 10"^ . . . 10^ Jy (dotted blue, bottom 
to top). Finally, the thick dotted black line is an empirical rela tion 
for the maximum SMBH binary mass versus redshift (Eq. IIOII . 




6F/€ (log Jy) 

Fig. 5. — Sources with orbit period modulation II: Signif- 
icant modulation: The number of point sources with flux modu- 
lation greater than the limiting sensitivity, versus that sensitivity. 
As FigureO except the number is shown versus SF/e, assuming a 
relative flux modulation 5 per orbit period given by Eq. 1151 

To this point we've considered the sensitivity limits for 
single measurements. By folding together multiple cycles 
we improve our sensitivity to faint sources. A priori, the 
flux sensitivity increases as the square root of the number 
of measurements Nobs', the figure of merit becomes 



SNR = (5 



F 



F„ 



(17) 
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Additionally, the survey duration Tabs and cadence AT = 
Tabs /Nobs naturally define a minimum and maximum 
range of periods to which the survey is sensitive: periods 
Pobs between approximately 2Tobs/Nobs and Tobs- This 
subtlety has minimal impact, as binaries cluster near the 
longest orbital periods. As implied by the ratio of Eqs. 
[T71 and simply rescaling the distribution of all bina- 
ries whose peak-to-trough modulations are above the de- 
tection threshold [Figure [5] produces the distribution of 
sources with detectable periodic flux modulations [Figure 

In the above we have assumed the flux is purely sinu- 
soidally modulated, suitable to weakly Doppler-boosted 
synchrotron emission detected in the radio. By contrast, 
a highly relativistic outflow will generally be strongly 
Doppler boosted, tightly beamed, and emit most power 
preferentially at high energies. As noted earlier, beamed 
emission from the rare jets pointing towards us could 
be detected farther away [Eq. [T3]. Few binaries have 
bright jets and precession periods less than a year [Fig- 
ure [2]. However, any transverse structure oc S6 in the 
narrow jet beam modulates the emitted flux on much 
shorter timescales Tmod — Sd/flp ~ y/ fbeamTprec- The 
modulation in high-energy flux with time {S{t)) will gen- 
erally contain a broad spectrum of frequencies, with fre- 
quency content depending on the details of the jet beam 
shape and its orientation relative to the line of sight. Be- 
cause many orders of magnitude more binaries have much 
longer precession periods, these modulations in principle 
allow surveys of short cadence to identify sources with 
long natural periods. To order of magnitude and focusing 
only on bright jets in the local universe, the total number 
of sources increases as (solid angle of jet) x (inspiral time 
starting at a 1 yr/59 precession period) x {d{l3) / 69)^ , the 
increased volume to which a beamed jet can be seen. 
Fixing /3 in the center and flnal factors by requiring 
Tprec = lyr/Se oc [Eq. (O] and substituting for d{l3), 
in principle the number of detectable precessing sources 
increases relative to the small-statistics tail Figure [2] by 
a substantial factor 

Nprec,beam ^ ea— 3/2 ^ f^3/4 
Tvr beam 

A much less favorable scaling applies at larger distances: 

following Eq. [TTl we expect N oc 60~^^'^ oc Jbelm- 
practice, however, the long-timescale nonperiodic modu- 
lations that should arise from slowly precessing, beamed 
AGN jets should be difficult to distinguish from back- 
ground AGN variability, as described in § |4] below. 

3.1. Ultrarelativistic outflows 

To this point wc have implicitly assumed the radio flux 
to be weakly modulated on the orbit or precession pe- 
riod, as expected from prompt emission along a weakly 
relativistic orbit. By contrast, if these outflows are sim- 
ilar to conventional ultrarelativistic outflow models for 
AGN jets, then time of flight, generation, and reprocess- 
ing delays may modify any intrinsic structure, delaying 
and distorting the expected pattern from precession or 
the orbit. 

Slowly processing ultrarelativistic jets share many 
qualitative features with more familiar ultrarelativistic 
outflows like short GRBs. In both cases, emission can 
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Fig. 6. — Sources with orbit period modulation III: Stack- 
ing significant modulation: For the small seed (SE) model, the 
cumulative number of sources A''(> SF^/W^ / e) at any instant 
which could be recovered by a survey with 1/d (thick) or 7/ yr 
(dotted), assuming one year observing time and stacking all Nobs 
flux measurements into a composite statistic. Each curve is di- 
rectly comparable to the dark b lue cu rve in Figure [31 after shifting 
the flux scale of the latter by •/Nobs- 
he shock-driven, powered by interaction with the ambi- 
ent medium. At best, emission from this shock will be 
delayed from prompt emission by time-of-flight and de- 
celeration delays. Synchrotron radiation at low frequen- 
cies can be self-absorbed; this reprocessing delays and 
distorts any modulations imposed by the orbit. Finally, 
only observers over a range of angles 9j + l/T can see 
the jet (for 6j the jet opening angle, F the bulk Lorentz 
factor); it may become visible only at late times. Just 
from outflow dynamics, the radio power expected from a 
relativistic outflow should plausibly peak on timescales 
of weeks to months after prompt emission. 

Therefore, ultrarelativistic outflows if present can fur- 
ther bias us towards the longest timescales: sources that 
produce the shortest timescales will average away, leav- 
ing the multiple- year modulations expected among wide 
SMBH binaries^ 

4. AGN BACKGROUNDS AND TARGETED SEARCHES 

To this point we have assumed a blind survey, lim- 
ited only by background and detector noise. In practice 
additional sources of noise intrinsic and extrinsic to the 
binary system can severely complicate its detection or 
idcntiflcation as a truly periodic source. On the contrary, 
highly relativistic SMBH binaries have unique kinematic 
and gravitational wave signatures that help distinguish 
true binary modulation from a background. Below we 
discuss two particular sources of noise (AGN variabilty 
and scintillation), and two methods to verify a partic- 
ular modulation has binary origin (multiple timescales 
and pulsar timing). 

4.1. AGN backgrounds and targeted searches 

The host systems of binaries might be identified 
by residual accretion-powered emission. For example, 
accretion- fueled AGN jets persist ©(Myr) after the 

By contrast, an extended epoch of delayed emission will signif- 
icantly enhance our prospects of detecting single brief events, such 
flares accompanying the merger of SMBH mergers; sec Kamble et 
al 2011. 
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black holes stop accreting and separate from the disk. 
Magnetically-powered jets may also already be identi- 
fied. Conceivably, a targeted search could much more 
efficiently identify variability candidates. 

On the other hand, AGN jets arc well-known to be 
strongly variable, with flux roughly pcrforraing a ran- 
dom walk on long timescales; see, e.g., iKrolikl (IT998I) and 
references therein. Likewise, these magnetically-driven 
flows could be variable. A periodic modulation super- 
imposed on a variable jet would be significantly more 
difficult to detect. 

To illuminate the targeting problem, we parameterize 
the worst possible scenario: every merging binary has a 
bright residual radio jet, emitting at a significant frac- 
tion of the binary's eddington luminosity {eoidLedd, with 
(old ^ e) and just as variable as active AGN. For time 
changes less than a y ear, AGN fluxes ar c roughly random 
walks in amplitude (jHughes et al.l [19921) Specifically, 
if w is a white random variable, a prototypical AGN flux 
and flux power spectrum is 



FAGN{t)- 



edd 



wdt 



SfU)^ / dte'^'f'' {FAGNiO)FAGNit)) 



( Coldic 



dd 



(/yr) ^yr 



(18) 



(19) 



normalizing to &crit relative variations on 1 year 
timescales. In this worst case hypothesis, for any de- 
tectable jet, AGN variability is the dominant noise. 
As both the signal and dominant noise propagate from 
the source to our detector, the flgure of merit for de- 
tectable sources depends only on ratios of properties at 
the sourceQ 

SNR=^= (20) 
^Sf(\IP)IT 



eq^{P/P^a.f{l~q^)/{l + q^) /Txlyr 



^oldSc 



p2 



eoldScrit 



X (2^AfG/(c/3„,a.)') '^'f , 

1 + q"^ 



(21) 



where for clarity we have replaced {3^ by the period [Eq. 

If indeed the residual jet radio power is brighter than 
modulation in the magnetically-driven radio emission 
(e < ^oid^crit), the factors in this a priori expression are 
potentially unfavorable to targeted searches. As black 
hole light crossing times [MG/c") are less than a few 
minutes, targeted searches for periodic AGN variability 
from a jet superimposed on an existing source would only 
be sensitive to short timescales (i.e., P < vTyrSOOs ~ 
1 d). Even then, detection could only occur if the binary 

We assume no correlation between AGN variability and the 
binary mass. 

Unlike the previous case, the AGN provides a unique noise 
realization for all detectors; repeated measurements at the same or 
multiple detectors does not improve our sensitivity. 



jet and residual jet fluctuations are finely tuned to nearly 
the same flux (i.e., e ~ eoidS)- 

On the other hand, this Poynting-dominatcd outflow 
isn't driven by accretion and is therefore likely far more 
stable than accretion-driven jets. Additionally, any resid- 
ual jets should have dimmed and expanded in the thou- 
sands to million years since accretion terminated. As- 
suming the transverse crossing time sets the characteris- 
tic AGN jet variability timescale, a not-face-on residual 
jet of age tj should have variability timescales > Ojtj-. 
as the jet ages, it loses the capacity to undergo rapid 
changes. As detailed treatment of this kind of residual 
AGN jet background is beyond the scope of this paper, 
we can neither endorse nor rule out the possibility that 
periodic variations might be accessible on top of likely 
jet variability backgrounds. 



4.2. Scintillation 

Emission from distant binaries propagates through 
plasma in their host, our galaxy, and the intergalac- 
tic medium. The motion of this intervening medium 
also introduces significant relative amplitude fiuctuations 
("scintillation"). Scintillation reprocesses both our tar- 
get modulated signal and any colocated background. 
The amount of scintillation depends strongly on the line 
of sight to the source, the plasma content and motion 
near the merging binary, the observing frequency, the 
angular size of t he sour c e, and t he modu l ation timescale 
of interest; see iWalked (|1998[ ). iRickettI (|1990[) and ref- 
erences therein. Because of the wide range in possible 
parameters that are unrelated to the intrinsic nature of 
the source, we do not model scintillation-induced back- 
grounds in detail here. 

To sonic degree, modeling is unnecessary. To the ex- 
tent that background AGN activity contaminates a sig- 
nal, a phenomenological approach to AGN activity au- 
tomatically includes all noise sources, scintillation in- 
cluded. Likewise, scintillation has characteristic power- 
law scal i ngs o f amplitude and timescale with frequency 
(|Walkeii |1998() . allowing identification and partial sub- 
traction of these modulations. 

That said, scintillation naturally produces 0(1) fiux 
variations on an unknown but long timescale comparable 
to any periodicity of interest. For extended sources ex- 
periencing strong refractive scattering, scintillation nat- 
urally occurs on long timescales, roughly the time an 
intervening plasma element v would need to cross the 
source size R. While the size of the emission region is 
poorly known, the lower bound of the horizon size insures 
this crossing time is large 

t.cin oc R/v ~ 2Oday(M/lO^Af0)(i?/Af)(5Okm/s/i;) 

where we scale v to typical galactic plasma velocities. In 
this case, typical scintillation flux variations should be 
large, with SF/F oc (^ ^0)^'^^^° ^ 0{1); based on Figure 
1 from iWalked ()1998[ ) we expect a transition frequency 
i/o between 1 to ten times a fiducial 1 GHz observing 
frequency, depending on the line of sight. For binary 
candidates found by radio surveys, multiband foUowup 
observations will be required to disentangle the impact 
of scintillation from any intrinsic modulation. 
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Fig. 7. — Typical lifetime: The average log lifetime (logT) ver- 
sus scaled flux F/e, including sources with observed orbital period 
< 1 yr. The brightest sources have short typical lifetimes, as they 
consist of the most massive nearby SMBH binaries in tight orbits. 

4.3. Unique time signatures versus confusing variability 
backgrounds 

As exemplified by OJ287 (|Katzl I1997D . an ob- 
served periodicity could be fit by other processes 
produced by single black holes, such a s a warped 
disk causing the black hole to precess (|Kat3 [Tail) 
or periodic accretion fr om or excited by a com- 
panion (see, e.g., 11.11 lArtvmowicz fc Lubowl [19961: 
IHavasaki et "S] l2007inFarris et al.l l20nr and references 
therein). Binary black holes can likewise produce pe- 
riodic emission through other mechanisms, including or- 
biting accretion flows and modes excited in the circumbi- 
nary disk. A detailed treatment of all possible back- 
grounds is beyond the scope of this paper. However, the 
mere presence of an unmodeled periodicity alone does 
not uniquely determine this mechanism. 

On the other hand, extremely sensitive surveys that re- 
cover many merging events may find a few with distinc- 
tive signatures: either (a) a distinctive ultrarelativistic 
(i.e., Dopplcr-boostcd and shapiro-dclay-distortcd) light 
curve; (b) multiple periodic timescales; or (c) a noticablc 
chirp. The most distinctive signature, a noticable chirp, 
occurs surprisingly often among detectable binaries. For 
example, an all-sky survey with F/e = O.lJy can po- 
tentially observe between 300 — 1000 binaries with less 
than 1 yr observed periods [Fig. [3] . That population has 
a ~ 10% probability of one member undergoing a chirp 
from 1 yr to 1 min within 10 yr. Moreover, the most mas- 
sive and brightest members generally have shorter chirp 
times; sec, e.g., the average log lifetime in Figure [7] 

4.4. Pulsar timing GW signature 

The massive, wide binaries which produce the most nu- 
merous and brightest jets are precisely the same sources 
whose gravitatio nal wave emission is targeted by pu l- 
sar timing arrays (jSesana et al.ll2009l : lHobbs et al.l(2010( ). 
Pulsar timing arrays are very likely to have a significantly 
shorter range than direct electromagnetic searches. To 
order of magnitude, the ratio of radio to gravitational 
wave power is ~ M'^r'^oj^ / Ledd, or 



L 



GW 



Lj rripC/ 



(27r/P)2/32 



On the other hand, the relative energy flux sensitivity 
of a gravitational wave survey with characteristic strain 
sensitivity hs 



F„. 



(23) 



As a result, the relative luminosity distance sensitivity 
of pulsar timing arrays and EM surveys goes as 



J'L.gw 
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(24) 
(25) 



3 X 10 




(j^.ad»o/mJy)(A//lO«A/0) ^^, 
(e/0.002) 



:3 X l(fe-^{P/yif{M/l(fMQ)f3'^ (22) 



where we ha ve scaled hf, to existi ng pulsar timing array 
sensitivities (jYardlev et al.l [20Tol ). Future pulsar tim- 
ing arrays will have enough reach to identify sources 
up to their confusion limit. Conservatively assum- 
ing only one sourc e per frequency can be distinguished 
(jSesana et al.ll2009l ) , at least a few sources over the entire 
sky could be individually resolved. Depending on how 
well the network localizes sources on the s ky, the confu- 
sion limit could be significantly lower (cf. iBovle fc PenI 
l2010f l. Nonetheless, barring extremely inefficient conver- 
sion of jet power to radio, direct electromagnetic sur- 
veys will have substantially greater sensitivity than pul- 
sar timing. A pulsar timing survey by itself will only 
be more effective if the efficiency e and flux F imply 
iV(> F/e) is small - in other words, only if radio surveys 
can find at best a handful of sources over the entire sky. 

Nonetheless, electromagnetic and gravitational wave 
surveys naturally complement each other. Owing to 
considerable uncertainties in the emission geometry and 
model, only a gravitational wave signature can con- 
firm a jet candidate corresponds to an SMBH binary. 
Conversely, starting with a well-localized binary candi- 
date from pulsar timing (jSesana fc Vecchiol[20im . elec- 
tromagnetic observation s can tightly localize a suitably- 
modulated companion (jSesana et al.l 1201 IbD with suit- 
able overall circumbinary spe ctra (jTanaka et al.l 120111 : 
IMilosavlievic fc Phinnevl[2005D . Subsequently, as obser- 
vations provide a candidate light curve corresponding to 
a known binary, archival and directed EM surveys can 
find more, at higher confidence. Eventually, in the limit 
of perfect coupling, a joint EM and gravitational wave 
survey could construct joint confidence in a proposed re- 
construction of its data as a superposition of nearby and 
distant binaries plus background variability. 

5. DISCUSSION 

Each SMBH binary should possess a pair of 
magnetically-powered bright jets attached to each hole 
during the last stages of gravitational-wave driven in- 
spiral. Supcrmassivc black holes therefore ubiquitously 
produce long-lived, modulated, and moderately bright 
emission even in (local) vacuum, without direct accre- 
tion. Though each jet is faint and though mergers are 
rare, gravitational-wave-driven inspiral is so slow that 
many modulation c ycles exis t per binary. Unlike rare 
merger events (jKaplan ct all [20ll[ ). these modulations 
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should easily be seen with fu ture blind radio surveys like 
VAST (jMurphv et al.|[20Tll ). only excepting pessimistic 
choices for the efficiencies eradio: eRriH and obscurat ion 
/geo- For example, following iPalenzuela et alJ ()2010bf ) in 
adopting the conservative value CEdd — 0.002, we predict 
30-60 observable binaries with observable orbital periods 
< lyr could be seen by the VAST survey (lO^deg sky 
coverage, and either 1/day or 7/year cadence; see Figure 
[6l). Alternatively , if we include the diffuse emission found 
bv lMoesta et al.l (|2011[ ). then eEdd is significantly larger, 
and roughly ten times more sources could be detected. 
Our results are insensitive to the precise merger rate we 
assume. 

At any time, the most accessible sources are the 
largest black hole binaries in the widest (~ year) orbits. 
Nonetheless, surveys with shorter cadence can achieve 
comparable sensitivity by stacking cycles. The gravita- 
tional wave emission fro m these systems is a lso the target 
of pulsar timing arrays (jSesana et al.ll2009() . Unless very 
little jet power escapes in modulated radio jets, radio sur- 
veys will generally be more efficient than pulsar timing 
surveys for the same binaries. 

The region surrounding a SMBH offers many oppor- 
tunities to produce variability on comparable timescales 
and comparable bands. In special but not exceptionally 
rare circumstances, this mechanism can be identified by 
multiple timescales - depending on the scale, some com- 
bination of effects from its ultrarelativistic orbit (Doppler 
and Shapiro), spin-orbit precession, and in exceptional 
circumstances the gravitational wave chirp. 

If most merging binaries accreted and produced long- 
lived AGN jets before decoupling from the circumbinary 
disk, then this long-lived remnant likely has significant 
variability which could severely limit the cffiacy of a ra- 
dio survey. Targeted searches towards known jets must 
overcome similar variability challenges. 

While we emphasize and scale fiducial results to radio- 
band measurements; our expressions involve only flux 
sensitivity limits and efficiencies and scale to any bandF^ 
On the other hand, our preferred model of sinusoidal 
modulations visible equally well in all directions likely 
scales to low-energy emission. High energy emission nat- 
urally arises from an ultrarelativistic and tightly beamed 
outflow, with dominant modulations from precession 
and subdominant modulations from orbital motion (e.g., 
from curvature near the source, for marginally-detected 
lines of sight into tight orbits). In that case, the proba- 
bility of detecting prompt emission from any one jet will 
be reduced by /Att. On the other hand, the internal 
structure of ultrarelativistic jets versus angle could mag- 
nify the effect of small changes in angle. A handful of 
tightly beamed jets could exhibit a distinctive rise and 
fall on the longest observing timescales. 

In this paper we only address the prospects for detect- 
ing a periodically modulated outflow in the radio, after 
the binary decouples completely from its circumbinary 
disk. We assume no accretion. However, so long as grav- 

At sufficiently high energies, loss than one photon arrives per 
orbit period in a reasonable-scale detector area. This effect is sig- 
nificant only for extremely low mass binaries in ultrarelativistic 
orbits, exceptionally high photon energies, or exceptionally low ef- 
ficiencies: M < O.5xlo3M0(d/Gpc)(E.^/keV)i/2.i;3/2/gi/2^ E^g^ 

for e ~ 10~® , the most interesting sources would still produce pho- 
ton modulations in the X-ray. 



itational radiation dominates the inspiral, our calcula- 
tions can be easily generalized to any partially- or totally- 
accretion-powered EM signal, just by replacing Eq. [T] 
for Lj{f3) and by modeling the modu lation d{P) in Eq. 
[Ml For example, iChang et all (|2010f ) propose an inner 
circumbinary disk will be accreted onto a more massive 
companion, with disk luminosity L oc T~5/4 ^ pio pg^k- 
ing at L ~ O.lLEdd- Observations will easily distinguish 
between these generalizations: each generically predicts 
different scalings of number versus flu x, period, and mo d- 
ulation shape. For example, Haiman et al.l (j2009al ) dis- 
cuss how the flux distribution N{> F) could be used to 
extract information about the relative scaling of inspiral 
and emission with orbit period. 

In our paper, we assume emission tracks the underlying 
SMBH dynamics; the outflow provides enough power for 
potentially significant reach in a quiet band (e.g., e > 
10~^); and reprocessing or relativistic kinematics at most 
weakly distorts the signal. The details of the emission 
mechanism remain unknown and merit further detailed 
investigations. 

For convenience, rather than discuss how the disk gen- 
erates a magnetic field and how the outflow produces 
prompt radio power, we have adopted a universal effi- 
ciency e = eEddCradio for all binaries, regardless of compo- 
nent masses, separation, or disk structure. Each product 
in the composite efficiency likely depends on binary and 
disk properties. For example, the efficiency eEdd is equiv- 
alent to the strength of the field threading the horizon; 
as described in Footnote [U this field strength depends 
on the disk structure and accretion flow. Likewise, the 
radio efficiency Cradio depends on the emission process, in 
turn possibly dependent on the amount of coronal mat- 
ter in the interior. Further investigations are needed to 
quantify the magnitude, variability, and parameter de- 
pendence of this composite efficiency. 

Finally, we have implicitly made strong assumptions 
about magnetic field order, particularly for the widest 
unequal-mass ratio binaries. The passage of the black 
hole through the magnetic field liberates the local instan- 
taneous energy density, possibly beamed along the in- 
stantaneous ambient field direction. So long as the mag- 
netic field is steady in magnitude and direction during all 
observations Tots, only the orbit's kinematics and out- 
fiow variability lead to modulations in the observed flux. 
However, for the widest orbits close to the circumbinary 
disk, each black hole will move through a spatially- and 
temporally- fluctuating field, causing additional variabil- 
ity. Detectable sources with such wide orbits should be 
rare [Figure[T]. Future MHD simulations of circumbinary 
disks should extract the correlation functions of magnetic 
fiux in the circumbinary region, to better determine the 
level of variability expected. 

DLK and AK are partially supported by NSF award 
AST-1008353. ROS is supported by NSF award PHY- 
0970074, the Bradley Program Fellowship (BCS), and 
the UWM Research Growth Initiative. AS is supported 
by the Max Planck Society. The authors thank M. Er- 
acleous, J. Krolik, H. Bignall, L. RezzoUa, P. Moesta, 
D. Alic, O. Zanotti, and J. Lazio for helpful feedback, 
as well as the anonymous referee for thoughtful sugges- 
tions. ROS and AS also thank Aspen Center for Physics, 



11 



where this work was completed, and all participants of its 2011 program on Galaxy and Supermassive Black Hole 

Coevolution. 



REFERENCES 



Apostolatos, T. A., Cutler, C, Sussman, G. J., & Thorne, K. S. 
1994, Phys. Rev. D, 49, 6274 

[ADS] 

Artymowicz, P. & Lubow, S. H. 1996, ApJ, 467, L77+ 
[ADS] 

Arun, K. C, Babak, S., Berti, E., Cornish, N., Cutler, C, Gair, 
J., Hughes, S. A., Iyer, B. R., Lang, R. N., Mandel, I., Porter, 
E. K., Sathyaprakash, B. S., Sinha, S., Sintes, A. M., Trias, M., 
Van Den Broeck, C., & Volonteri, M. 2009, Classical and 
Quantum Gravity, 26, 094027 
[ADS] 

Beckwith, K., Hawley, J. F., & Krolik, J. H. 2009, ApJ, 707, 428 

[ADS] 

Bertone, S. & Conselice, C. J. 2009, MNRAS, 396, 2345 
[ADS] 

Bloom, J. S. & et al. 2009, arXiv:0902.1527 

[ADS] 

Bode, T., Bogdanovic, T., Haas, R., Healy, J., Laguna, P., & 

Shoemaker, D. 2011, arXiv:1101.4684 [URL] 
Bogdanovic, T., Eracleous, M., & Sigurdsson, S. 2009, ApJ, 697, 

288 

[ADS] 

Bogdanovic, T., Reynolds, C. S., & Miller, M. C. 2007, ApJ, 661, 
L147 

[ADS] 

Boroson. T. A. & Lauer, T. R. 2009, Nature, 458. 53 

[ADS] 

Boyle, L. Sz Pen, U.-L. 2010, ArXiv e-prints 
[ADS] 

Burke-Spolaor, S. 2010, MNRAS, 1574 
[ADS] 

Chang, P., Strubbe, L. E., Menou, K., & Quataert, E. 2010, 
MNRAS, 407, 2007 
[ADS] 

Chen, X., Madau, P., Sesana, A., & Liu, F. K. 2009, ApJ, 697, 
L149 
I [ADS] I 

Comerford, J. M., Gerke, B. F., Newman, J. A., Davis, M., Yan, 
R., Cooper, M. C, Faber, S. M., Koo, D. C, Coil, A. L., 
Rosari o, D. J., & Dutton, A. A. 2009, ApJ, 698, 956 
flADSll 

Demorest, P., Lazio, J., Lommen, A., Archibald, A., 
Arzoumanian, Z., Backer, D., Cordes, J., Demorest, P., 
Ferdman, R., Freire, P., Gonzalez, M., Jenet, R., Kaspi, V., 
Kondratiev, V., Lazio, J., Lommen, A., Lorimer, D., Lynch, R., 
McLaughlin, M., Nice, D., Ransom, S., Shannon, R., Siemens, 
X., Stairs, L, Stinebring, D., Reitze, D., Shoemaker, D., 
Whitcomb, S., & Weiss, R. 2009, in Astronomy, Vol. 2010, 
AGB Stars and Related Phenomenastro2010: The Astronomy 
and Astrophysics Decadal Survey, 64 — h 

Dotti, M., Montuori, C, Decarli, R., Volonteri, M., Colpi, M., & 
Haardt, F. 2009, MNRAS, 398, L73 
[ADS] 

Eracleous, M., Boroson, T. A., Halpern, J. P., & Liu, J. 2011, 

(arXiv:1106.2952) [URL] 
Farris, B. D., Liu, Y. T., & Shapiro, S. L. 2011, Phys. Rev. D, 84, 

024024 

[ADS] 

Haiman, Z., Kocsis, B., & Menou, K. 2009a, (arXiv;0904.1383) 
[URL] 

Haiman, Z., Kocsis, B., Menou, K., Lippai, Z., & Frci, Z. 2009b, 
Classical and Quantum Gravity, 26, 094032 
[ADS] 

H ayasaki . K., Mineshige, S., & Sudou, H. 2007, PASJ, 59, 427 



Hobbs, G., Archibald, A., Arzoumanian, Z., Backer, D., Bailes, 
M., Bhat, N. D. R., Burgay, M., Burke-Spolaor, S., Champion, 
D., Cognard, I., Coles, W., Cordes, J., Demorest, P., Desvignes, 
G., Ferdman, R. D., Finn, L., Freire, P., Gonzalez, M., Hessels, 
J., Hotan, A., Janssen, G., Jenet, F., Jessner, A., Jordan, C, 
Kaspi, v., Kramer, M., Kondratiev, V., Lazio, J., Lazaridis, K., 
Lee, K. J., Levin, Y., Lommen, A., Lorimer, D., Lynch, R., 
Lyne, A., Manchester, R., McLaughlin, M., Nice, D., Oslowski, 
S., Pilia, M., Possenti, A., Purver, M., Ransom, S., Reynolds, 
J., Sanidas, S., Sarkissian, J., Sesana, A., Shannon, R., 
Siemens, X., Stairs, I., Stappers, B., Stinebring, D., Theureau, 
G., van Haasteren, R., van Straten, W., Verbiest, J. P. W., 
Yardley, D. R. B., & You, X. P. 2010, Classical and Quantum 
Gravit y, 27, 084013 
[Ml 

H ughes, P. A., AUer, H. D., & Aller. M. F. 1992, ApJ, 396, 469 

H ughes, S. A. 2002, MNRAS, 331, 805 
flADST 

Jenet, 1'., Finn, L. S., Lazio, J., Lommen, A., McLaughlin, M., 
Stairs, I., Stinebring, D., Verbiest, J., Archibald, A., 
Arzoumanian, Z., Backer, D., Cordes, J., Demorest, P., 
Ferdman, R., Freire, P., Gonzalez, M., Kaspi, V., Kondratiev, 
v., Lorimer, D., Lynch, R., Nice, D., Ransom, S., Shannon, R., 
& Siemens, X. 2009, ArXiv e-prints 
I [ADS] 

Johnston, S. et al. 2007, PASA, 24, 174 
[ADS] 

Kaplan, D. L., O'Shaughnessy, R., Sesana, A., & Volonteri, M. 
2011, A pJ, 734, L37-I- 
'[ADS]' 

Katz, J. I. 1997, ApJ, 478, 527 
[ADS] 

King, A. R., Pringle, J. E., & Livio. M. 2007, MNRAS, 376, 1740 
[ADS] 

Komossa, S. 2006, Mem. Soc. Astron. Italiana, 77, 733 

[ADS] 

Komossa, S., Burwitz, V., Hasinger, G., Predehl, P., Kaastra, 
J. S., & Ikebe, Y. 2003, ApJ, 582, L15 
[ADS] 

Krolik, J. 1998, Active Galactic Nuclei (Princeton Series in 

Astrophysics) 
Lawrence, A. & Elvis, M. 2010, ApJ. 714. 561 

[ADS] 

Loeb, A. & Gaudi, B. S. 2003, ApJ, 588, L117 
[ADS] 

Lyutikov, M. 2011, Phys. Rev. D, 83, 064001 
pDS]J 

Marscher , A. P. & Jorstad, S. G. 2011. ApJ. 729. 26 
'[ADS]' 

Megevand, M., Anderson, M., Frank, J., Hirschmann, E. W., 
Lehner, L., Liebling, S. L., Motl, P. M., & Neilsen, D. 2009, 
Phys^Rev. D, 80, 024012 
'pLDS] 

MilosavljevicS, M. & Phinney, E. S. 2005, ApJ. 622. L93 

[ADS[ 

Moesta, P., Alic, D., Rczzolla, L., & Zanotti, O. 2011, arXiv 
preprints 

Murphy, T., Chatterjee, S., et al. 2011, PASA, in prep. 

Ofek, E. O., Frail, D. A., Breslauer, B., Kulkarni, S. R., Chandra, 

P., Gal- Yam, A., Kasliwal, M. M., & Gehrels, N. 2011, ApJ, 

arXiv:1103.3010 

|[ADS] 

O'NeiU, S. M., Miller, M. C, Bogdanovic, T., Reynolds, C. S., & 
Schnittman, J. D. 2009, ApJ, 700, 859 

r[ADS] 

Palenzuela, C, Lehner, L., & Liebling, S. L. 2010a, Science, 329, 
927 
[[ADS] 

— . 2010b , Science, 329, 927 



12 



Pcssah, M. E., Chan, C, & Psaltis, D. 2006, Physical Review 
Letters. 97, 221103 
[ADS] 

Rickctt, B. J. 1990. ARA&A, 28, 561 
[ADS] 

Rodriguez, C, Taylor, G. B., Zavala, R. T., Peck, A. B., Pollack, 
L. K., & Romani, R. W. 2006a, ApJ, 646, 49 

[ADS] 

— . 2006b, ApJ, 646, 49 
[ADS] 

Roedig. C, Dotti, M., Sesana, A., Cuadra, J., & Colpi, M. 2011, 
MNRAS, 979 
[[ADS] 

Rossi, E. M., Lodato, G., Armitage, P. J., Pringle, J. E., & King, 
A. R. 2010, MNRAS, 401, 2021 
[ADS] 

Rothstein. D. M. & Lovelace, R. V. E. 2008, ApJ, 677, 1221 
[ADS] 

Schnittman, J. D. 2010, ArXiv e-prints 

[ADS] 

Schnittman, J. D. & Krolik, J. H. 2008, ApJ, 684, 835 |[URL]| 
Sesana, A., Gair, J., Berti, E., & Volonteri, M. 2011a, 
Phys. Rev. D, 83. 044036 

[ADS] 

Sesana, A., Roedig, C., Reynolds, M. T., & Dotti, M. 2011b, 
ArXiv e-prints 

[ADS] 

Sesana, A. & Vecchio, A. 2010, Phys. Rev. D, 81, 104008 

[ADS] 

Sesana, A., Vecchio, A., & Volonteri, M. 2009, MNRAS, 394, 2255 
[ADS] 

Sesana, A., Volonteri, M., & Haardt, F. 2007, MNRAS, 377, 1711 
[ADS] 

Shen, Y. & Loeb, A. 2010, ApJ, 725, 249 

[ADS] 

Sillanpaa, A., Haarala, S., Valtonen, M. J., Sundelius, B., & 
Byrd, G. G. 1988, ApJ, 325, 628 



Smith, K. L., Shields, G. A., Bonning, E. W., McMuUen, C. C., 
Rosari o, D. J., & Salviander, S. 2010, ApJ, 716, 866 
I [ADS] I 

Stone, N. & Loeb, A. 2011, MNRAS, 224 
[ADS] 

Tanaka, T. & Menou, K. 2010, ApJ. 714, 404 

[ADS] 

Tanaka, T., Menou, K., & Haiman, Z. 2011, ArXiv e-prints 
[ADS] 

Tsalmantza, P., Decarli, R., Dotti, M., & Hogg, D. W. 2011, 
ArXiv e-prints 

U rry, C M. & Padovani, P. 1995, PASP, 107, 803 
lADST 

Valtonen, M. J., Lehto, H. J., Nilsson, K., Heidt, J., Takalo, 
L. O., Sillanpaa, A., Villforth, C., Kidger, M., Poyner, G., 
Pursimo, T., Zola, S., Wu, J.-H., Zhou, X., Sadakane, K., 
Drozdz, M., Kozicl, D., Marchev, D., Ogloza, W., Porowski, C., 
Siwak, M., Stachowski, G., Winiarski, M., Hentunen, V.-P., 
Nissinen, M., Liakos, A., & Dogru, S. 2008, Nature, 452, 851 
IIADS] 

Volonteri, M. 2010, A&A Rev., 18, 279 

Volonteri, M. & Begelman, M. C. 2010, MNRAS, 409, 1022 
[ADS] 

Volonteri, M., Haardt, F., & Madau, P. 2003, ApJ, 582, 559 
[ADS] 

Volonteri, M., Lodato, G., & Natarajan, P. 2008, MNRAS, 383, 
1079 

[ADS] 

Walker, M. A. 1998. MNRAS, 294, 307 
[ADS] 

Wegg, G. & Bode, J. N. 2010, arXiv;1011.5874l[URL]| 
Yardley, D. R. B., Hobbs, G. B., Jenet, F. A., Vcrbicst, J. P. W., 
Wen, Z. L., Manchester, R. N., Coles, W. A., van Straten, W., 
Bailes, M., Bhat, N. D. R., Burke-Spolaor, S., Champion, D. J., 
Hotan, A. W., & Sarkissian, J. M. 2010, MNRAS, 407, 669 
[Ml 



